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Abstract 



Recent measurements from the HERMES and SMC collaborations show a 
remarkably large azimuthal single-spin asymmetries Ajjl and Ajjt of the proton 
in semi-inclusive pion leptoproduction j*{q)p — > ttX. We show that final-state 
interactions from gluon exchange between the outgoing quark and the target 
spectator system lead to single-spin asymmetries in deep inelastic lepton-proton 
scattering at leading twist in perturbative QCD; i.e., the rescattering correc- 
tions are not power-law suppressed at large photon virtuality at fixed Xbj. 
The existence of such single-spin asymmetries requires a phase difference be- 
tween two amplitudes coupling the proton target with = it^ to the same 
final-state, the same amplitudes which are necessary to produce a nonzero pro- 
ton anomalous magnetic moment. We show that the exchange of gauge particles 
between the outgoing quark and the proton spectators produces a Coulomb- 
like complex phase which depends on the angular momentum of the pro- 
ton's constituents and is thus distinct for different proton spin amplitudes. The 
single-spin asymmetry which arises from such final-state interactions does not 
factorize into a product of distribution function and fragmentation function, 
and it is not related to the transversity distribution Sq{x,Q) which correlates 
transversely polarized quarks with the spin of the transversely polarized target 
nucleon. 



1 Introduction 

Single-spin asymmetries in hadronic reactions have been among the most difficult 
phenomena to understand from basic principles in QCD. The problem has become 
more acute because of the observations by the HERMES [Q] and SMC collabora- 
tions of a strong correlation between the target proton spin Sp and the plane of the 
produced pion and virtual photon in semi-inclusive deep inelastic lepton scattering 
ip^ — s> i'nX at photon virtuality as large as = 6 GeV^. Large azimuthal single- 
spin asymmetries have also been seen in hadronic reactions such as pp^ nX 0, 
where the target antiproton is polarized normal to the pion production plane, and in 
pp A^X 0], where the hyperon is polarized normal to the production plane. 

In the target rest frame, single-spin correlations correspond to the T-odd triple 
product iSp-p^^ x q, where the phase i is required by time-reversal invariance. The dif- 
ferential cross section thus has an azimuthal asymmetry proportional to IjvI |^sin6'q^sin(/) 
where is the angle between the plane containing the photon and pion and the plane 
containing the photon and proton polarization vector S'p. In a general frame, the az- 
imuthal asymmetry has the invariant form -^e^varP^SpP'^^q^ where the polarization 
four- vector of the proton satisfies 5^ = — 1 and S'p ■ P = 0. 

In order to produce a correlation involving a transversely-polarized proton, there 
are two necessary conditions: (1) There must be two proton spin amplitudes 
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M[7*p(t/p) F] with Jp = ±i which couple to the same final-state \F >; and (2) 
The two amplitudes must have different, complex phases. The correlation is propor- 
tional to Im(M[Jp^ = +i]*M[JJ = -i]). The analysis of single-spin asymmetries thus 
requires an understanding of QCD at the amplitude level, well beyond the standard 
treatment of hard inclusive reactions based on the factorization of distribution func- 
tions and fragmentation functions. Since we need the interference of two amplitudes 
which have different proton spin = ±| but couple to the same final-state, the or- 
bital angular momentum of the two proton wavefunctions must differ by AL^ = 1. The 
anomalous magnetic moment for the proton is also proportional to the interference 
of amplitudes M['y*p{Jp) — > F] with = ±| which couple to the same final-state 
\F >. 

Final-state interactions (FSI) in gauge theory can affect deep inelastic scattering 
reactions in a profound way, as has been demonstrated recently |0 . The rescattering 
of the outgoing quark leads to a leading twist contribution to the deep inelastic cross 
section from diffractive channels 7*p — > qqp', and the interference effects induced 
by these diffractive channels cause nuclear shadowing. Here we shall show that FSI 
also provide the required phase needed to produce single-spin asymmetries in deep 
inelastic scattering. 

The dynamics of the constituents in the target can be described by its light-front 
wavefunctions, ipn/pi^i, k±i, Xi), the projections of the hadronic eigenstate on the free 
color-singlet Fock state |n > at a given light-cone time r = t + z/c. The wavefunctions 
are Lorentz-invariant functions of the relative coordinates Xi = kf / = {k^ + 

kf)/{P^+P^) and k±i [with J2i=i = I and J27=i = O-i]' they are independent 
of the bound state's physical momentum P"*" and P± 0. The physical transverse 

momenta are p±i = XiP± + k^i. The Aj label the light-front spin projections of the 
quarks and gluons along the quantization z direction. If a target is stable, its light- 
front wavefunction must be real. Thus the only source of a nonzero complex phase in 
leptoproduction in the light-front frame are final-state interactions. The rescattering 
corrections from final-state exchange of gauge particles produce Coulomb-like complex 
phases which, however, depend on the proton spin. Thus M[Yp{Jp = ±|) ^ F] = 
\M[^*p{Jp = ±i) — > F] \ e*^±. Each of the phases is infrared divergent; however the 
difference Ax = x+ ~ X- is infrared finite and nonzero. The resulting single-spin 
asymmetry is then proportional to sinA^. 

2 A Model Calculation of Single-Spin Asymme- 
tries in Gauge Theory 

We shall calculate the single-spin asymmetry in semi-inclusive electroproduction 7*p — 
HX induced by final-state interactions in a model of a spin-| proton of mass M with 
charged spin-^and spin-0 constituents of mass m and A, respectively, as in the QCD- 
motivated quark-diquark model of a nucleon. The basic electroproduction reaction is 
then 7*p — > q{qq)o, as illustrated in Figs. |I| and ^ We shall take the case where the 
detected particle H is identical to the quark. One can take the asymmetry for a de- 
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Figure 1: The final-state interaction in the semi- inclusive deep inelastic lepton scat- 
tering tp^ — > £VX. 



tected hadron by convoluting the jet asymmetry result with a realistic fragmentation 
function; e.g. Dq^T^x^z^Q"^). 

The amplitude for the 7*p (l{(l(l)o can be computed from the tree and one- 
loop graphs illustrated in Fig. A spin asymmetry will arise from the final-state 
interactions of the outgoing charged lines. The = +| two-particle Fock state is 
given by 0, |] 



^ two particle (-^ ) ^1 

d^fc I dx 



0, 



x{l — x)167r'^ 



; xP+ 

2 ' 



+-: xP+, ki 



where 



^l,{x,k^) = {M+f)^, 

The scalar part of the wavefunction (p depends on the dynamics. In the perturbative 
theory it is simply 



ip = ip{x, k±) 



n-x 



M2 



(3) 



In general one normalizes the Fock state to unit probability. 
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Figure 2: The tree (a) and one-loop (b) graphs for 7*p (l{(l<l)o- The interference of 
the two amphtudes with = ±1/2 provides the proton's single-spin asymmetry. 

Similarly, the = two-particle Fock state has components 

7' - (4) 

2 

The spin-flip amplitudes in (Q) and have orbital angular momentum projection 
= +1 and —1 respectively. The numerator structure of the wavefunctions is char- 
acteristic of the orbital angular momentum, and holds for both perturbative and 
non-perturbative couplings. 

We require the interference between the tree amplitude of Fig. ^ and the one loop 
graph of Fig. 0d. The contributing amplitudes for 7*p have the following 

structure through one loop order: 

{M+'^)C{h + ^'-f^g,) (5) 
C(ft + .fl,.) (6) 
-) C {h + t^g2) (7) 





-T) 




-T) 




H) 







' A ' ' Stt 

{M + ^)C{h + z'-^g,), (8) 
lA on 



where 



C = - c/ ei P+ 2 A (1 - A) (9) 

1 

n +A(l-A)(-M2 + i^ + - 



h = rr, 2 T2- ■ (10) 



The quark light-cone fraction A = is equal to the Bjorken variable Xbj up to 
corrections of order 1/Q. The label it / JJ- corresponds to = ±|. The second label 
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Figure 3: The light-cone frame used is p = {p~^,p~,'p±) = (P+, M^/P+, 0_l) and 
q = {q~^,q~,q±) with q^ = 0. The virtual photon and produced hadron define the 
production plane which we will take as the z — x plane. 



t / i gives the spin projection = ±| of the spin-i constituent. Here ei and 62 
are the electric charges of q and (gg)o, respectively, and g is the coupling constant of 
the proton-g-(gg)o vertex. The first term in (|^) to is the Born contribution of the 
tree graph. The crucial result will be the fact that the contributions gi and (72 from 
the one-loop diagram Fig. ^ are different, and that their difference is infrared finite. 
A gauge particle mass Xg will be used as an infrared regulator in the calculation of 
gi and g2. 

The calculation will be done using light-cone time-ordered perturbation theory, or 
equivalently, by integrating Feynman loop diagrams over dk~ . The light-cone frame 
used is p = {p~^,p~',p±) = (P"*", M^/P"*", 0_l) and q = {q^,q^,q±) with g+ = and 
q~ = 2q ■ p/P~^, q± = Qx with = —q^. The Bjorken variable is A = Q'^/2q ■ p = 
Q'^/2Mu. Since q~^ = 0, light-cone time-orderings where the virtual photon produces 
a qq pair do not appear. 

The light-cone formalism is invariant under boosts in the z direction: P^ 7P+. 
It reduces to a laboratory frame when P"*" = M. If we take q to lie in the z—x plane in 
this frame, q = (g^, g^, g^) = (Q, 0, — z/); i.e., q is oriented at an angle 9iab = tan~^ ^ 

from the negative z direction. This is illustrated in Fig. S. Here i' is the laboratory 
energy of the photon. In the Bjorken scaling limit with and u large, and A = Xbj 
fixed, the angle 6iab — > 0, so the light-cone laboratory frame and usual laboratory 
frame with q taken in the —z direction are identical. 

The covariant expression for the four one-loop amplitudes of diagram Fig. [pD is: 



^one-loop (11) 



ig ej 62 



X 



- m2 + ie) {{k + qf - m? + te){{k - r)^ - + ie){{k - Pf - + it) 



2(27r)4 J X X {x- A) (1 - x) 
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X 



1 

(A2+(A;x-rx)2)-i. 



{k- 



P- 



where we used k^ 



{x~A)P+ J \^ ' (l~x)P+ 

xP^ . The numerators N'{I) are given by 



2P+\/A X (M +—) q- 

X 

2P+VA X {- ) q- 



2P+VA X 



-k^-ik'^ 



X 



2P+yA X {M+—) q 



X 



(12) 

(13) 

(14) 
(15) 



where q = = . For the [current]- [gauge propagator]- [current] factor, in 

Feynman gauge only the —g^~ term of the gauge propagator —g'^'^ contributes in the 
Bjorken hmit, and it provides a factor proportional to q~ in the numerator which 
cancels the q~ in the denominator of the gauge propagator. Therefore the result 
scales in the Bjorken limit. 

The integration over k~ in (0) does not give zero only if < x < 1. We first 
consider the region A < x < 1. 

^one— loop 



(16) 



P+dx 



P+4 XX (x-A) {1 



X) 



X 



X 



^ {l-x)P+ xP+ / 1 (l-x)P+ 



(m.2+(fcx+g±)^)-ie 
xP+ 



P- 



{l-x)P+ 



(A2 + (fcx-rx)^)-»£ 
{x-A)P+ 



The result is identical to that obtained from light-cone time-ordered perturbation 
theory. 

The phases Xi needed for single-spin asymmetries come from the imaginary part 
of (0), which arises from the potentially real intermediate state allowed before the 
rescattering. The imaginary part of the propagator (light-cone energy denominator) 
gives 



-in 6 \ P 



(A^ + kl) 
(1 -x)P+ 

1 A2 



(m^ + {k^_ + gl)2) 



xP+ 



^ ^(^ - ^ - ^) 
r qj_ 



(17) 



where 

5 = 2A^^1%^. (18) 
Q± 

Since the exchanged momentum 5P~^ is small, the light-cone energy denominator 

corresponding to the gauge propagator is dominated by the (x-a)~^^ term. This 

gets multiplied by (x — A), so only {k± — f±Y + \g appears in the propagator, in- 
dependent of whether the photon is absorbed or emitted. The contribution from the 
region < x < A thus compliments the contribution from the region A < x < 1. 

We can integrate ([1^) over the transverse momentum using a Feynman parame- 
trization to obtain the one-loop terms in (|^) to (|]). 

/•I 1 

Qi = da — 2 r-? — (19) 

a(l-a)rl + aA2 + (l-a)A(l- A)(-M2 + ^ + ^) ^ ' 







1 a 



Jo a(l - a)rl + aX^ + (1 - a)A(l - A)(-M2 + !^ + ^) ' 

Although not necessary for our analysis, we will assume for convenience that the final- 
state interactions generate a phase when exponentiated, as in the Coulomb phase 
analysis of QED. The rescattering phases e*-^' {i = 1,2) with Xi = ^^^"^(^sf^^) 
thus distinct for the spin-parallel and spin-antiparallel amplitudes. The difference in 
phase arises from the orbital angular momentum k± factor in the spin-flip amplitude, 
which after integration gives the extra factor of the Feynman parameter a in the 
numerator of (72- Notice that the phases Xi are each infrared divergent for zero gauge 
boson mass Ag ^ 0, as is characteristic of Coulomb phases. However, the difference 
Xi ~ X2 which contributes to the single-spin asymmetry is infrared finite. We have 
verified that the Feynman gauge result is also obtained in the light cone gauge using 
the principal value prescription. The small numerator coupling of the light-cone gauge 
particle is compensated by the small value for the exchanged = 5P^ momentum. 

The virtual photon and produced hadron define the production plane which we 
will take as the z — x plane. The azimuthal single-spin asymmetry transverse to the 
production plane is given by 



e,e2 2(AM + m)ri ^2 ^2 

^ ^ fl+A(l-A)(-M=' + — + 



AM + m ) + 



A 1-A' 

L 1 j.y± I 1 1 L \ I I I 

rl A(l-A)(-Af2 + f 

The linear factor of = reflects the fact that the single spin asjnnmetry is pro- 
portional to ■ g X r where q ~ —vz and Sp = ±y. Here A = x^j. 

Our analysis can be generalized to the corresponding calculation in QCD. The 
final-state interaction from gluon exchange has the strength — > (^^^^^(/i^). The 
scale of a, in the MS scheme can be identified with the momentum transfer carried 
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by the gluon /i^ = e'^^^'^{k±_ — rj_)^ [Q. The matrix elements of the proton to its 
constituents will have the same numerator structure as the perturbative model since 
they are determined by orbital angular momentum constraints. The strengths of the 
proton matrix elements can be normalized by the anomalous magnetic moment and 
the total charge. In QCD, is the magnitude of the momentum of the current quark 
jet relative to the virtual photon direction. Notice that for large r_L, Vy decreases as 
as{r^)xi,jMrx_ inr^ ^ physical protou mass M appears since it is present in the ratio 

of the Lz = 1 and = matrix elements. This form is expected to be essentially 
universal. 



3 Model Predictions 

We show the predictions of our model in Fig. ^ for the asymmetry Vy = A^J^'^ of 

the iSp ■ q X Pq correlation based on Eq. (PT|). As representative parameters we 
take as = 0.3, M = 0.94 GeV for the proton mass, m = 0.3 GeV for the fermion 
constituent and A = 0.8 GeV for the spin-0 spectator. The single-spin asymmetry 
Vy is shown as a function of A at = 0.5 GeV in Fig. 4a and as a function of r± 
at A = 0.15 in Fig. 4b. The Hermes asymmetry A^^jf contains a kinematic factor 
K = ^y/l — y = \f^^\J^^ because the proton is polarized along the direction of 
the incident electron. The resulting predictions for KVy are shown in Figs. 4c and 
4d. Note that f = pq — q is the momentum of the current quark jet relative to the 
photon direction. The asymmetry as a function of the pion momentum requires a 
convolution with the quark fragmentation function. 



4 Summary 

We have calculated the single-spin asymmetry in semi-inclusive electroproduction 
induced by final-state interactions. We have shown that the final-state interac- 
tions from gluon exchange between the outgoing quark and the target spectator 
system leads to single-spin asymmetries in deep inelastic lepton-proton scattering 
at leading twist in perturbative QCD; i.e., the rescattering corrections are not power- 
law suppressed at large photon virtuahty at fixed xi,j. The azimuthal single- 
spin asymmetry Vy transverse to the photon-to-pion production plane decreases as 
as{r\)xhjMr i_\\.\ir ^/ for large rj_, where r±_ is the magnitude of the momentum of 
the current quark jet relative to the virtual photon direction. The fall-off in instead 
of compensates for the dimension of the q-q -gluon correlation. The mass M of 
the physical proton mass appears here since it determines the ratio of the = 1 and 
= matrix elements. We have estimated the scale of as as 0{r\). The nominal 
size of the spin asymmetry is thus CFOis{r\)ap where is the proton anomalous 
magnetic moment. 

It is usually assumed that the cross section for semi-inclusive deep inelastic scat- 
tering at large factorizes as the product of quark distributions times quark frag- 



mentation functions [10, ITT|. Our analysis shows that the single-spin asymmetry 
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(b) 





Figure 4: Model predictions for the single spin asymmetry of the proton in electro- 
production resulting from gluon exchange in the final state as a function of A = Xi,j 
and quark transverse momentum r±. The parameters are given in the text. 



which arises from final-state interactions does not factorize in this way since the re- 
sult depends on the < plip^Aiplp > proton correlator, not the usual quark distribution 

derived from < p\ipq{^)il!q{0)\p > evaluated at equal light-cone time = 0. In par- 
ticular, the spin asymmetry is not related to the transversity distribution Sq{x, Q) 
which correlates transversely polarized quarks with the spin of the transversely po- 
larized target nucleon. 

Our results are directly applicable to the azimuthal correlation of the proton spin 
with the virtual photon to current quark jet plane, which can be deduced from jet 
measures such as the thrust distribution. The sin0 correlation of the proton spin 
with the photon-to-pion production plane as measured in the HERMES and SMC 
experiments can then be obtained using the usual fragmentation function. Detailed 
comparisons with experiment will be presented elsewhere. Our approach can also be 
applied to single-spin asymmetries in more general hadronic hard inclusive reactions 
such as e+e~ — > h}X and pp — > h}X. 
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